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It has recently been pointed out that global solutions of Einstein's equations for a 3-brane universe 
embedded in 4 spatial dimensions give rise to a Friedmann equation of the form H (x p on the brane, 
instead of the usual H oc ^/p, which is inconsistent with cosmological observations. We remedy this 
problem by adding cosmological constants to the brane and the bulk, as in the recent scenario of 
Randall and Sundrum. Our observation allows for normal expansion during nucleosynthesis, but 
faster than normal expansion in the very early universe, which could be helpful for electroweak 
baryogenesis, for example. 
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During the past year, much has been written about 
the possibility of having compactified extra dimensions 
with large radii [Q. In the original proposal, Mp was 
related to the radius 60 of the N compact dimensions 
by Mp = M^(M6o)^, where M is the new fundamen- 
tal quantum gravity scale, which could in principle be 
as low as 1 TeV. If so, this would be a partial solution 
of the hierarchy problem, i.e., why the weak scale, Mw, 
is 17 orders of magnitude smaller than the Planck scale, 
Mp: it is because bo is, for some reason, much larger 
than M~^. If 60 3> Mp^, as is necessary if M ~ M\y, 
the particles and fields of the standard model must be 
restricted to stay on a 3-dimensional slice (brane) of the 
full iV-|-3 spatial dimensions; otherwise particle propaga- 
tion in the new dimensions would already have been seen 
in accelerator experiments. But even with the restriction 
of the brane, the idea implies many possibly observable 
effects at accelerators. It also poses severe challenges for 
cosmology. In this letter we will address one of the cos- 
mological problems, and comment upon an unexpected 
connection to the question of precisely how the hierarchy 
problem is solved using the extra dimensions. 

Our starting point is the observation recently made 
by Binetruy, Deffayet and Langlois [||^ that the Fried- 
mann equation for the Hubble expansion rate of our 3D 
universe is modified, even at very low temperatures, by 
the presence of an extra dimension, y, compactified on 
a circle or an orbifold. Allowing for the possibility of a 
cosmological constant Aft in the full 4 spatial dimensions, 
called the bulk, the new Friedmann equation for the scale 
factor a of our brane is fGUil 
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instead of the usual relation, H = ^ ptjiMp. pt is the 
total (vacuum plus matter) energy density on the brane. 
This expression is derived, as will be explained below, 
from the 5D action 



(2) 



S= I d^x dy ^\ {iM'^n - Aft + £,,_,) 
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where the action, >Cbrano, for the matter living on the 
brane results in a stress-energy tensor parametrized as 
Tj^ = J(6y)diag(-/9t,pt,pt,pt,0). An interesting aspect 
of this result is that fact that, in order to find consistent 
global solutions to the Einstein equations in the (4-1- 1)- 
D spacetime, it is necessary to add a second brane 0], a 
mirror of our own, having equal and opposite energy den- 
sity. This topology can be motivated from string theory. 
In the Hof ava-Witten picture [g[ of the nonperturbative 
regime of the E^ x E^ string theory, the string coupling is 
interpreted as an eleventh compact dimension with a Z2 
symmetry that truncates the spectrum in order to keep 
only sixteen supercharges in lOD, i.e., an A/" = 1 super- 
symmetry in 4D after compactification on a Calabi-Yau 
manifold. There is good evidence that over a wide 
range of energies the theory behaves like a 5D theory 
compactified on a Z2 orbifold with two 3-branes, viewed 
as the remanants of the lOD hypersurfaces where the E^ 
gauge groups were living. The two 3-branes can also be 
seen as D3-branes of the type / string theory ||I| . 

Naively, one would expect that at distances much big- 
ger than the size of the fifth dimension, the effects of 
the extra compact dimension become small corrections 
to the usual 4D equations; thus when H~^ 3> 6o> one 
should recover the standard cosmology. However the 
presence of the mirror brane contradicts this logic. Let 
us choose the range of the compact coordinate to be 
ye [—1/2, +1/2]. The solutions of the Einstein equations 
for the scale factor a{y) behave ||g| like ao(l -I- A|j/|/2), 
with A ^ Pt. Because the points y = ±1/2 are iden- 
tified, the derivative is discontinuous at this point, and 
a" /a = A{S{y) — 5{y — 1/2)). The Einstein equations 
identify the delta functions with the energy densities of 
the two respective branes. In the limit as &o ~* 0, the two 
branes overlap, and their energy densities cancel to first 
order in pt because they are equal and opposite. There- 
fore only terms of order {a' /a)'^ ^ A^ ^ p^ survive, even 
at arbitrarily late times in cosmological history. The re- 
sulting expansion rate (1) is probably incompatible with 
big bang nucleosynthesis, which is extremely sensitive to 
how the Hubble rate varies with the the energy density. 



hence temperature. Even if one tunes M so that the 
altered expansion rate (1) still gives the correct helium, 
it is likely that the other elements will come out wrong, 
since their rates of production depend quite differently 
on the temperature. 
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Figure 1. Qualitative dependence of the 3D scale factor a{y) 
on the compact dimension y in the solutions of; (solid line) ref. 
[0], with vanishing bulk cosmological constant; and (dashed 
line) ref. H, with At, given by eq. (H). a(0) is nonzero but 
exponentially small in the latter. 

From eq. (Q), one can imagine a very simple escape 
from this dilemma 0. Suppose there is a cosmological 
constant A localized on our brane (and correspondingly 
—A on the mirror brane, although this value will be cor- 
rected by terms of order p in the presence of matter on 
the branes), so that pt — K + p, where p now denotes 
the energy density of normal matter or radiation on the 
brane, as opposed to vacuum energy. One can choose Af, 
to exactly cancel the A^ terms in eq. (|lj) , and furthermore 
fix the value of A in terms of M and Adp : 
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where ± refers to the two respective branes. Condition 
(3) insures the cancellation of the effect of A;, by A^ in (1), 
whereas (4) adjusts the overall rate of expansion to agree 
with the usual result. The new Friedmann equation then 
becomes the conventional one, plus a correction which is 
quadratic in the density: 
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We have distinguished the values of p on the two branes 
by the subscript to emphasize that they need not~in fact, 
cannot-be the same. The brane with the positive solution 
has a rate of expansion that is consistent with all current 
cosmological observations as long as the normal rate has 
been recovered by the epoch of nucleosynthesis, which 



will be true if < p+ < 0.1 (1 MeV)-^ 
finds the constraint that 

M > 10 TeV, 



< A. One thus 



(6) 



which is not much more severe than other accelerator and 
astrophysical limits that have recently been placed on the 
new gravity scale. The other brane must have p_ < 0, 
since otherwise H^ < 0, which has no solution. 

The condition (||) is precisely what is needed to get 
a static universe in the case of vanishing p: the nega- 
tive cosmological constant in the bulk cancels the posi- 
tive A^ from either brane. The solutions to the Einstein 
equations in this case were recently studied by Randall 
and Sundrum (RS) |^, but for very different reasons: 
they found that the weak scale hierarchy problem is nat- 
urally solved on one of the branes, even if M '^ Mp, and 
6o ^ 50Mp^. This comes about because the metric ten- 
sor has an exponential dependence on the coordinate of 
the compact 5th dimension (see Figure 1). Using the line 
element 

ds^ = -n'^{t,y)dt^ + a^{t,y)Sijdx'dx^ +b{t,yfdy^, (7) 

it is straightforward to verify the time-independent solu- 
tion 
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One then observes that, even if all mass parameters in 
the Lagrangians for matter on the branes are of the or- 
der Mp, the physical masses on the brane at y = 1/2 
are suppressed by the factor e"*^", which can be of order 
Mw /Mp with only a moderate hierarchy between bo and 
Mp^ ~ M~^. Since g^i, enters differently in the kinetic 
than the mass terms for a scalar field, once the kinetic 
terms are canonically normalized, masses get multiplied 
by a(l/2) ^ e~^l'^ . This idea therefore appears to be 
a much more natural solution to the hierarchy problem 
than the original proposal, which required h^M to be of 
order (J^'lpjM)'^!^ , where N is the number of extradi- 
mensions. 

We now see that the static solution of RS is the start- 
ing point for our idea, which is to recover the normal 
expansion of the 3D universe by perturbing large, bal- 
ancing cosmological constants in the bulk and the branes 
by a small density of matter or radiation on the branes. 
Intuitively, it is clear that solutions with nonvanishing p 
must exist, but we will now take some time to demon- 
strate this explicitly, in the vicinity of our brane. We 
were not able to find global solutions in closed form once 
matter with an arbitrary equation of state p = ujp was 
introduced. However, we are really most interested in 
the expansion rate on our own brane, so it suffices to 
solve the Einstein equations in that region. To simplify 
the appearance of the solutions, we will translate the y 



coordinate by i/ ^ j/ + 1/2, so that the brane which we 
inhabit is located at y = 0. 

We must solve the Einstein equations for the metric 
(^, now allowing for time dependence in a, 6 and n. It is 
always possible to chose a gauge so that n{t, 0) is constant 
at y = 0, without introducing g^^ elements in the metric. 
We will make this choice, and drop all terms involving h 
since they are not relevant for the solution in the imme- 
diate vicinity of the brane. With this simplification, the 
5D Einstein equations, G^^ = M~'^T^y, become 
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in the vicinity of y = 0. Close to our brane, the nonzero 
elements of the 5D stress-energy tensor are 

Too = n^{p + I^)S{by) + n^At + V{b)) ; 
Tu - a^p - A)5{hy) - a^Af, + V{h)) ; 
T55^-h\At + V{b) + V'{h)/h) (14) 

where S(by) = b~^S{y) is the generally covariant form of 
the delta function. There are also source terms propor- 
tional to h~^d{y — 1/2) at the mirror brane, but these 
will not directly concern us in what follows. The terms 
involving V{h) would result if there is a potential that 
stabilizes the compact dimension. Their presence does 
not qualitatively change any of our conclusions, so we 
set V{h) to zero in what follows. 

The generalization of the static solution (0) can be 
parametrized as 



a{t, y) = ao(i) exp(i^|j/| -I- \A2y^ 
6(t,y)=6oexp(iB|y| + iB22/' 
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n(t,y)=exp(iAA|y|-KiAA2j/2 
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By our choice of gauge for time, there is no no(t) function. 
We have not assumed separability of the solution here, 
since the coefficients A,B,J\f, etc., need not be static; 
however we will see that their time-dependence arises en- 
tirely from that of p and p. The fact that 60 is constant 
in time is not obvious, but will be proven to be consistent 
with eqs. (^fsj). 

As in ref. gVthe linear-in-|y| coefficients, A and A/", are 
determined by the singular parts of eqs. (||) and (|ll]), i.e., 



those involving the delta functions and second spatial 
derivatives. One finds that 
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Therefore, to obtain solutions that are growing in the di- 
rection of the mirror brane, as are needed to solve the 
hierarchy problem on our own, we would have to chose 
A < here, about which we shall say more below. The 
analogous coefficient B is not determined in this way be- 
cause h" appears nowhere in the Einstein equations. But 
it is constrained by eq. (p!3|). Inserting the ansatz ( [l5| ) in 
this equation, and taking lj — pj p to be constant (which 
is a weak restriction since p and p refer only to the matter 
and radiation), one can eventually show that 
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and that it is consistent to take 60 = 0. Thus the scale 
factor of the compact dimension, although it expands 
inside the bulk, is strictly constant on our brane. Eq. 
(O) is not a complete specification for B since A2 and 
Ml are not yet known, but in fact we will never need B 
for determining the Friedmann equation on our brane. 

It remains to satisfy the nonsingular parts of the other 
Einstein equations, (g - |l^), near y = 0. The knowledge 
of A and M is all that is needed to specify eq. (|l^) at 
y = because no second derivatives appear. One obtains 
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where the second equation follows from using our previ- 
ous determination of A and Af,, eqs. (§-§)• The leading 
term reproduces the usual prediction of general relativ- 
ity, and the second term corresponds to the quadratic 
correction in eq. (ra). Indeed, in light of the energy con- 
servation law on the brane, p = —3H(p + p), which is 
true regardless of the extra dimension |g] , (g) is the only 
relation consistent with ( |l8[ ) when u! — p/ p is assumed 
to be constant. 

In contrast to the new Einstein equation (|2|) associ- 
ated with the 5th dimension, the Goo and Gu equations 
(^ and (^1]) depend on the quadratic coefficients A2 and 
M2 at y = 0, because of the presence of a" and n" . With 
two equations in two unknowns, it is always possible to 
find values of A2 and ^2 such that the resulting equations 
for ao(t) are consistent with (p[) and (Eq). Therefore eqs. 
(0) and (O) add no new information on the brane, al- 
though they would be necessary if one wanted to deduce 
the full y dependence of the solutions in the bulk. 

In the above derivation, it was shown that the brane 
whose masses are small by the RS mechanism must have 



A < 0. Unfortunately, we already saw in eq. (|5|) that 
the brane with negative A must have an energy density 
P- < 0, which is not the case in our universe. This would 
appear to be a serious problem for the RS idea. However, 
see the "Note added," below. 

The conditions &m for the brane and bulk cosmolog- 
ical constants look strange at first, so some words of mo- 
tivation are in order. Although when p = 0, A-|_ — — A_ 
on the two branes, when p ^ 0, a global solution to the 
Einstein equations is needed in order to derive the exact 
relation between A_|_ and A_ when p 7^ 0, since it involves 
all the coefficients of the expansion ( p^ ) pO[ | : 
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In addition to this topological relation derived from the 
spacetime geometry, there is also a relation involving At,. 
We argue that the latter is a stringent consistency con- 
dition similar to the global tadpole cancellation in string 
theory (see for instance ||ll[]): 
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i.e. the global effective cosmological constant must van- 
ish. In the solution of RS, this condition reduces to 
A§ + QM^Ki, = 0, which is the relation they needed to 
obtain a global solution to Einstein equations. 

The condition (|2^) can be understood if the cosmo- 
logical constants are viewed as an effective description of 
the Ramond-Ramond fields of the underlying string the- 
ory: for instance the value of the (j) -t- l)-form to which 
a p-brane is coupled is reinterpreted as a cosmological 
constant on the p-brane. The condition ( [20[ ) will now be 
necessary to cancel the UV divergences of the string the- 
ory. The connection between the phenomenological sce- 
nario of RS and string theory has recently been examined 
by Verlinde [12| and his analysis concludes that the expo- 
nential dependence of the metric in the compact direction 
is identified with the renormalization group scale when 
using the AdS/CFT correspondence, which also corrob- 
orates the stringy origin of the RS mechanism. 

Since the normal expansion rate of the universe is 
only known to have held between nucleosynthesis and 
the present epoch (as was stressed in reference |13J) it 
would be interesting if the quadratic corrections to the 
new Friedmann equation (pi) started to become important 
above temperatures of several MeV. In the most natural 
version of the RS scenario, M and Mp are of the same 
order, so the corrections become important only at the 
Planck scale. However it is still a logical possibility to 
imagine that the fundamental scale M is much smaller 
than Mp. In this case one recovers the Arkani-Hamed 
et al. result that Mp = M^bo, which combined with 
the gravitational tests that restrict 60 ^ 1 mm, gives the 
constraint M > 10^ GeV. With such a large value of 



M, departure from normal expansion occurs only above 
temperatures T >, 1 TeV, which is not far above the elec- 
troweak scale. An intriguing possibility would be to in- 
crease the rate of expansion during the electroweak phase 
transition. If this occurred, standard model sphaleron in- 
teractions could easily be out of equilibrium in the broken 
phase |l3[, making electroweak baryogenesis more feasi- 
ble. 
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Note added: After acceptance of this work, we discov- 
ered ||l^ that the problem of the wrong-sign expansion 
rate at the second brane can be solved if the extra dimen- 
sion is taken to be noncompact, as suggested by ref. [[l6[ . 
Ref. [|5| shows that by considering multiple intersecting 
branes, the whole construction can be extended to any 
number of extra dimensions. 
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